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This letter reports new results from the MINOS experiment based on a two-year exposure to 
muon neutrinos from the Fermilab NuMI beam. Our data are consistent with quantum mechanical 
oscillations of neutrino flavor with mass splitting | Am 2 1 = (2.43 ±0.13) x 10" 3 eV 2 (68% confidence 
level) and mixing angle sin 2 (29) > 0.90 (90% confidence level). Our data disfavor two alternative 
explanations for the disappearance of neutrinos in flight, namely neutrino decays into lighter particles 
and quantum decoherence of neutrinos, at the 3.7 and 5.7 standard deviation levels, respectively. 

PACS numbers: 14.60.Lm, 14.60.Pq, 29.27.-a, 29.30.-h 



Several experiments [l[ have produced compelling evi- 
dence of the disappearance of neutrinos of a given lepton 
flavor. In previous publications Q, the MINOS exper- 
iment has also presented evidence for energy-dependent 
disappearance of muon neutrinos produced by the NuMI 
facility at Fermilab. Based on the number of events, that 
result provides evidence of the disappearance of at a 
level of 5.2 standard deviations. Such observations sup- 
port the description of neutrinos via two independent 
basis sets, the mass and the flavor eigenstates, with the 
bases related by the 3x3 PMNS matrix 0]. They im- 
ply that at least two of the neutrino eigenstates have 
non-zero mass. In this letter we present results obtained 
from a larger dataset than that used in Q. These re- 
sults provide a precision measurement of the oscillation 
parameters and further disfavor two other theoretical in- 
terpretations of neutrino flavor disappearance 0, Q . 

The MINOS detectors [|| and the NuMI beam line Q 
are described elsewhere. In brief, NuMI is a conventional 
two-horn-focused neutrino beam with a 675 m long decay 
tunnel. The horn current and position of the hadron 
production target relative to the horns can be configured 
to produce different energy spectra. MINOS consists 
of two detectors: a 0.98 kt Near Detector (ND) 1.04 km 
from the NuMI target; and a 5.4 kt Far Detector (FD) 
735 km from the target. Both are segmented, magnetized 
calorimeters that permit particle tracking. The curvature 
of muons produced in z^ + Fe — » fj,~ + X interactions [l9[ 
is used for energy determination of muons that exit the 
detector and to distinguish the component of the beam 
from the 6% intrinsic 77^ contamination. The energy of 
muons contained in the detector is measured via their 
range. Oscillations of into other neutrino flavors result 
in an energy-dependent depiction of interactions in 
the FD relative to the expectation based upon the ND 
measurement. 

The present letter describes results from data recorded 
between May 2005, and July 2007. Over this period, a 
total of 3.36 x 10 20 protons on target (POT) were accu- 
mulated for this analysis. A 1.27 x 10 20 POT subset of 
this exposure (hereafter referred to as Run I) forms the 
data set from Ref Q ■ In Run I and for most of the new 
running period (Run II), the beam line was configured to 
enhance production with energies 1-5 GeV (the low- 



energy configuration). An exposure of 0.15 x 10 20 POT 
was accumulated with the beam line configured to en- 
hance the i/fj, energy spectrum at 5-10 GeV (the high- 
energy configuration). The Run II data were collected 
with a replacement target of identical construction due to 
failure of the motion system of the first target. The new 
target was found to be displaced longitudinally ~1 cm 
relative to the first target, resulting in a 30 MeV shift in 
the neutrino spectrum. This effect is incorporated in the 
Monte Carlo simulation, and the Run I and Run II data 
sets are analyzed separately to account for this shift. 

The simulation of neutrino production and detection 
is accomplished with a model of hadron production in 
the target using Fluka Q and a GEANT3 H simulation 
of the beamline and detector. NEUGEN3.5.5 [n|, tuned 
to data from previous bubble chamber neutrino exper- 
iments and experiments with pion beams scattering on 
iron, is used to model neutrino interactions. As in our 
previous analysis, the Monte Carlo (MC) simulation of 
the neutrino flux was constrained to agree with the neu- 
trino energy spectrum in the ND collected in nine differ- 
ent configurations of the NuMI beam , thereby reduc- 
ing uncertainties in the flux prediction at the FD. Fig.[T] 
compares the simulation to the ND data acquired in the 
two configurations used in the oscillation analysis. 

Neutrino interactions in the MINOS detectors can ei- 
ther be charged-current, + Fc — > /i~ + X , or neutral- 
current, + Fc — > Up + X. In this analysis, only the for- 
mer are used because they identify the interacting neu- 
trino flavor and because the reconstructed energy best 
measures the full neutrino energy. To select charged- 
current events, we have implemented a new algorithm (11 1 
based on a multivariate likelihood including four vari- 
ables that characterize a muon track: the event length; 
the average pulse height per plane along the track; the 
transverse energy deposition profile of the track; and the 
fluctuation of the energy deposited in scintillator strips 
along the track. The new selection algorithm, along with 
a new track-finding algorithm, improve our efficiency to 
identify and select charged-current interactions in the FD 
from 75.3% using the previous selection to 81.5% in 
the current selection, in the absence of oscillations. The 
new selection reduces the neutral-current contamination 
in the charged-current sample from 1.8% in our previous 
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FIG. 1: Energy spectra in the MINOS ND for two of the nine 
beam configurations before and after tuning the Monte Carlo 
simulation to the ND data. The data combine Run I + Run 
II. Both configurations are utilized in the oscillation analysis. 



publication to 0.6% in the present analysis. The present 
analysis uses a larger fiducial mass of 4.17 kt in the FD, 
an increase of 2.9% over the mass used in Q. 

The measured energy spectrum at the ND is used to 
predict the energy spectrum at the FD. As in our previ- 
ous analysis , we compute a transfer matrix to correct 
for ~20% differences expected in the shape of the energy 
spectrum in the FD relative to the ND that arise from 
meson decay kinematics and from beam-line geometry 
0, EH ■ We have cross-checked this technique by compar- 
ison to other calculations of the FD spectrum [2J . 

The FD energy spectra were inspected only after 
the analysis procedure was finalized and basic data in- 
tegrity checks were performed. We observe 848 events 
in the FD for all energies 0-120 GeV produced by the 
NuMI beam, compared to the unoscillatcd expectation of 
1065±60 (syst.). In the low energy configuration alone, 
the number of events observed in the data is 730, to be 
compared with an expectation of 936 ± 53 (syst.). The 
observed energy spectrum of the events from the low- 
and high-energy datasets is shown along with the pre- 
dicted spectrum in Fig. [2] and the ratio of these data to 
the expected spectrum is shown in Fig. [3] 

Under the assumption the observed deficit is due to 
— > v T oscillations [l[ , a fit is performed to extract the 
parameters Am 2 and sin 2 (29) [131 ] using the expression 



1 - sin 2 (26»)sin 2 1.27Am 



E 



(1) 



where L[km] is the distance from the target, £?[GeV] is 
the neutrino energy, and |Am 2 | is measured in eV 2 . The 
FD data from Run I, Run II, and the high-energy run are 
separately fit to Eq. [TJ The best-fit parameters minimize 
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FIG. 2: Comparison of the FD data (points, with statistical 
uncertainties) from the low- and high-energy configurations 
with the predictions for the energy spectra with and with- 
out the effect of oscillations. The estimated neutral current 
(NC) background is indicated. 



the x 2 expression given in Q. The predicted oscillated 
spectrum includes the contamination from v T produced 
in the oscillation process. 

The effects of systematic uncertainties were evaluated 
by fitting modified MC in place of data. TableQ] gives 
the differences between the fitted values obtained with 
the modified and unmodified MC. The largest effects are: 
(a) the ±10.3% uncertainty in the absolute hadronic en- 
ergy scale, which is the sum in quadrature of a ±5.7% 
error in the calorimeter response to hadrons as derived 
from test beam measurements 1J], a ±2.3% uncertainty 
in the energy scale calibration, and a ±8.2% uncertainty 
in the simulation of neutrino production of hadrons in 
iron nuclei; (b) the ±3.3% relative uncertainty in the 
hadronic energy scale between the ND and FD; (c) the 
±4.0% uncertainty on the predicted FD event rate which 
is the sum in quadrature of the uncertainties on the de- 
tectors' fiducial mass, event selection efficiency and the 
POT counting; (d) the ±50% uncertainty on the neu- 
tral current contamination in the charged current event 
sample; and (e) the uncertainty on the muon momenta 
measured via range (±2.0%) or curvature (±3.0%). 

In fitting the data to Eq.Q] sin 2 (20) was constrained 
to lie in the physical region. To reduce the effect of the 
dominant systematic uncertainties in Table[T]((a) and (c) 



for | Am | , and (d) for sin 2 (29)) these three systematic 
uncertainties were included as nuisance parameters in 
the fit. The resulting best fit to the neutrino energy 
spectrum is shown in Fig. [2] and Fig. [3] We obtain 
|Am 2 | =(2.43 ±0.13) x 10" 3 eV 2 and sin 2 (20) > 0.95 at 
68% confidence level (CL.) The fit x 2 =90 for 97 de- 
grees of freedom. The resulting 68% CL. (Ax 2 =2.30) 
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FIG. 3: Ratio of the FD data and the expected spectrum in 
the absence of oscillations. Also shown are the best fit curve 
to Eq.[T]and the best fit to alternative models of neutrino dis- 
appearance [3, [1| . For display purposes, the data have been 
rebinned and the estimated oscillated NC background is sub- 
tracted. 



Uncertainty 


|Am 2 | 


sin 2 (2t9) 




(KT 3 eV 2 ) 




(a) Abs hadronic E scale (± 10.3%) 


0.052 


0.004 


(b) Rel hadronic E scale (± 3.3%) 


0.027 


0.006 


(c) Normalization (± 4%) 


0.081 


0.001 


(d) NC contamination (± 50%) 


0.021 


0.016 


(e) fi momentum (range 2%, curv 3%) 


0.032 


0.003 


[l)o v {E v < 10 GeV) (±12%) 


0.006 


0.004 


(g) Beam flux 


0.010 


0.000 


Total Systematic Uncertainty 


0.108 


0.018 


Expected Statistical Uncertainty 


0.19 


0.09 



TABLE I: Sources of systematic uncertainties in the mea- 
surement of |Am 2 | and sin 2 (2#). The values are the average 
shifts for varying the parameters in both directions without 
imposing the sin 2 (2(9) < 1 constraint on the fit. Correlations 
between the systematic effects are not taken into account. 
The dominant uncertainties are incorporated as nuisance pa- 
rameters in the fit of our data to Eq. [1] so as to reduce their 
effect on the oscillation parameter measurement (see text). 



present data set. Our measurement is the most precise 
determination of the mass splitting |Ato 2 |. 

Fig. 2] also shows that the previous MINOS result Q 
is in good agreement with the current measurement. 
Taken alone, the Run II data yield | A?n 2 | = (2.32±°;^g) x 
l(T 3 eV 2 and sin 2 (2(9)= 1.0, to be compared with 
(2.57to 2 2 3 Q ) x 10~ 3 eV 2 and sin 2 (26»)=1.0 from Run I. The 
two results are consistent at 68% C.L. We note that the 
value of 2.57x 10 -3 eV 2 for Run I differs from that quoted 
in Q because of our improved reconstruction and selec- 
tion of charged-current events and improved MC simula- 
tion of neutrino interactions. 

We have also fit the FD energy spectra to alternative 
models that have been proposed to explain the disappear- 
ance of neutrinos in flight, namely, the decay of neutrinos 
to lighter particles (Eq. 13 of [ij]), and the decoherence of 
the neutrino's quantum-mechanical wave packet (Eq. 5 
of [5j). Fig. [3] shows the ratios of the energy spectra aris- 
ing from our best fits to these alternative models to the 
prediction of the FD spectrum in the absence of dis- 
appearance. The x 2 f° r the best fit to the decay model 
is 104/97 d.o.f., while that for the decoherence model 
is 123/97 d.o.f. Given the A X 2 = 14 and 33 of these 
two models relative to the oscillation hypothesis, these 
models are disfavored with respect to the oscillation hy- 
pothesis at the 3.7 and 5.7 standard-deviation levels. 

In summary, we have presented updated measure- 
ments of neutrino oscillation parameters from the MI- 
NOS experiment. Based upon an exposure of 3.36 x 
10 20 POT from the NuMI beam, we obtain |Am 2 | = 
(2.43 ± 0.13) x 10" 3 eV 2 (68% C.L.) and mixing angle 
sin 2 (26») > 0.90 (90% C.L.). As the dataset presented 
here includes the subset analyzed in 0, these results su- 
persede our previous publication. Our data disfavor two 
alternative explanations for disappearance of neutrinos 
in flight, namely neutrino decays [4| into lighter particles 
and quantum decoherence of neutrinos Q at the 3.7 and 
5.7 standard-deviation level, respectively. 

This work was supported by the US DOE; the UK 
STFC; the US NSF; the State and University of Min- 
nesota; the University of Athens, Greece; and Brazil's 
FAPESP and CNPq. We are grateful to the Min- 
nesota Department of Natural Resources, the crew of the 
Soudan Underground Laboratory, and the staff of Fermi- 
lab for their contribution to this effort. 



and 90% C.L. (Ax 2 =4.61) intervals for the oscillation pa- 
rameters | Am 2 1 and sin 2 (26*) are shown in Fig.[2][l6|]. The 
MC predicts negligible backgrounds of 0.7 events from 
cosmic ray muons, and, at the best-fit value for |Am 2 | 
and sin 2 (2t9), 2.3 events from neutrino interactions in the 
upstream rock, 5.9 neutral current and 1.5 v T events in 
the final sample. If the fit is not constrained to the phys- 
ical region, |Am 2 |=2.33 x 10" 3 eV 2 and sin 2 (26»)= 1.07, 
with a 0.6 unit decrease in \ 2 - Correspondingly, the con- 
tours in Fig. @] are smaller than those expected for the 
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